A Delay Stochastic Process with Applications to Molecular Dynamics

Molecular processes of cell differentiation often involve reactions with delays. We develop a mathematical model that provides a basis for a rigorous theoretical analysis of these processes as well as for direct simulation. A discrete, stochastic approach is adopted because several molecules appear in small numbers only. Our model is a non-Markovian stochastic process. The main theoretical results include a constructive proof of the existence of the process and a derivation of the rates for initiation and completion of reactions with delays. These results guarantee that the stochastic process is an accurate and consistent description of the molecular system. They also serve as a theoretical justification of recent work on delay stochastic simulation. We apply our model to an important process in developmental biology, the formation of somites in the vertebrate embryo. Simulation of the molecular oscillator controlling this process reveals major differences between stochastic and deterministic models.